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SUMMARY 


A Microwave Ice Accretion Measurement Instrument (MIAMI) 
incorporating a microprocessor has been developed that can 

1. Detect the presence of ice and sound an alarm. 

2. Measure and digitally display ice thickness. 

3. Measure and digitally display ice accretion 
rate . 

4. Plot ice thickness, accretion rate and other 
parameters on a pen recorder. 

5. Log and print out a permanent record of ice 
thickness and ice accretion rate versus time. 

6. Store data for future statistical analysis and 
printout . 

Ice thickness indications to approximately 1/4 inch were 
frequently observed; however, it is believed that calibration 
techniques were good only up to 1/8 inch. Theory, verified 
during this program, indicates that larger ice thickness measure- 
ments are possible by resorting to larger transducers and lower 
operating frequencies. Ice thickness to 1 inch is not considered 
unreasonable . 

The instrument is non-intruslve , that is, it does not use 
probes; but is mounted under the surface being monitored and does 
not interfere with its aerodynamic properties. 

One or more transducers may be mounted anywhere on an airfoil 
or other aircraft surfaces - all can be monitored by a single micro- 
processor . 

Electro-thermal deicing of the MIAMI to remove existing ice 
is feasible. Ice removal is not necessary, however, until the ice 
thickness exceeds the dynamic range of the instrument. 

The MIAMI transducer is constructed from solid state components 
all of which have very little mass and are therefore capable of 
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being miniaturized and ruggedized. The unminiaturized experi- 
mental model was only 7" x 4” x 1 /i" in size. 

The development described herein represents research that was 
performed over a period of one year at IDEAL RESEARCH INC. The 
fundamental principle exploited is that layers of ice accreting on 
the surface of a resonant surface waveguide will alter its reso- 
nant frequency. The degree of shift of the resonance is related to 
thickness of the ice layers. To take advantage of this phenomenon 
the research proceeded in two areas: (l) the microwave transducer 

and (2) the microprocessor for monitoring resonant frequency and 
displaying ice thickness and accretion rate. 

The phenomenon observed and utilized are the result of the 
molecular and dielectric properties of water and ice. The work 
thus represents, for the first time, research into one of the 
basic physical phenomenon underlying the formation and adhesion of 
ice to aircraft surfaces. 
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DEVELOPMENT^ AND TES T OF A MICROWAVE ICE ACCRETION 
MEASUREMENT INSTRUMENT (MIAMI). 

1.0 INTRODUCTION 

The purpose of this program was to perform an engineering 
study to confirm that it is technically and economically feasible 
to develop a microwave instrument system (MIAMI) for research 
purposes that can simultaneously perform the following functions. 

a. Detect the presence of ice. 

b. Measure ice thickness. 

c. Measure ice accretion rate. 

"Icing instrumentation (Reference 1) is an area of research 
that is so fundamental that it is probably the most important area. 
There can not be an "icing science," without the ability to mea- 
sure the icing environment accurately. In the 1940's and 1950's, 
relatively clumsy and slow measurement techniques were used with 
questionable accuracy. Advances in instrumentation have been so 
dramatic since that time that it is now possible to achieve mea- 
surements of icing parameters with real time data displays. The 
"Microwave Ice Accretion Measurement Instrument" (MIAMI) developed 
during this program is one such device that can be used as a re- 
search, as well as an operational, instrument. 


1. Milton A. Beheim "Executive Summary of 
Aircraft Icing Specialists Workshop" NASA 
Conference Publication 2086 (FAA-RD-78-109) 
Aircraft Icing, July 19-21, 1978, Page 8. 
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2.0 RESULTS OF EXPERIMENTAL PROGRAM 


2.1 Description of the MIAMI 

The MIAMI developed in this program is a device that is 
responsive to the thickness of ice accreting on its surface. 

An illustration of the MIAMI transducer mounted on an airfoil is 
shown in Figure 1. Its real time output, illustrating the growth 
of ice on an airfoil mounted in the icing research tunnel (IRT), 
is shown in Figure 2. A photo of the MIAMI experimental model is 
shown in Figure 3. It consists of two major elements, a trans- 
ducer for mounting under the surface being monitored and a micro- 
processor connected to the transducer by a long multiconductor 
cable. The microprocessor sends signals to the transducer and 
based upon the response received from the transducer computes the 
ice thickness and accretion rate which are digitally displayed on 
the front panel. The displays are instantaneous. The MIAMI can 
be programmed to log and printout a permanent record of ice thick- 
ness and accretion rate as a function of time. A sample printout 
appears as Figure 4. It can also be programmed to store ice thick- 
ness and accretion rate in computer memory for later statistical 
analysis for use in icing cloud studies. 

Additional important features of the MIAMI transducer are 
that it is non-intrusive ; it does not use probes protruding into 
the air stream but is flush mounted, with the surface being moni- 
tored. The MIAMI is constructed from components that can be 
packaged as miniature integrated circuits having almost negligible 
mass, enabling it to be mounted almost anywhere on the skin of the 
aircraft including rotor blades and engine inlets. 

The transducer mounted just below the surface being moni- 
tored is a resonant microwave device, the resonant frequency of 
which is very sensitive to the thickness of ice accreting on its 
surface. As the ice layer grows, the resonant frequency is re- 
duced as illustrated in Figure 5 a and b. 
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FIGURE 



1 - MIAMI MOUNTED ON A TWO-DIMENSIONAL AIRFOIL 
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FIGURE 2 - TYPICAL PEN RECORDING OF MIAMI OUTPUT 
ILLUSTRATING ICE GROWTH IN REAL TIME 
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FIGURE 4 


ItUM NUMBER I 9 

TUNNEL PRAMSi VEL-100 TEHP<F>-30 DROPLET SI2E*12 WATER CONTENT-1.3 
COMMENTS* 

TRANSIT CDEPFICXENTt ;,93 


TIME 

D-INDEX 

TRANSIT ICE 

MIAMI ICE 

34.80 

367. 00 

18.60 

14.59 

61.20 

646.00 

37.20 

35.29 

94.76 

759. 00 

55.80 

45.99 

127.36 

937. 00 

74.40 

65.27 

163. 16 

1049.00 

93. 00 

79.88 

190.92 

1165.00 

111.60 

99. 07 

9. 00 

0.00 

0. 00 

0. 00 


8^25^81 RUN NUMBER: 10 

TUNNEL PRAMS: VEL-100 TEMP<F>-30 DROPLET S12E-15 WATER CONTENT-1. 3 

comments: 


TRANSIT 

COEFFICIENT: 

.93 


TIME 

D- INDEX 

TRANSIT ICE 

MIAMI ICE 

77.64 

512.00 

13.60 

24. 17 

117.24 

736.00 

37.20 

43.72 

154. 00 

865. 00 

55.80 

57. 05 

187.64 

1016.00 

74.40 

75.27 

216.20 

0.00 

0. 00 

0.00 

0. 00 


8^25^81 RUN number: 11 

TUNNEL PRAMS: VEL-100 TEMP<F>»30 DROPLET SIZE-19 WATER CONTENT-2.5 

comments: 

TRANSIT coefficient: .93 


TIME 

D-INLEX 

TRANSIT ICE 

MIAMI ICE 

36.92 

332. 00' 

18.60 

12.70 

87.80 

624. 00 

37.20 

33-34 

122.48 

745. 00 

55. SO 

44.61 

146.96 

846. 00 

74.40 

54.98 

172-52 

869. 00 

93. 00 

57.49 

191.92 

954. 00 

111.60 

67.32 

214.28 

928. 00' 

130.20 

64.20 

0. 00 

0.00 

0.00 

0. 00 


- SAMPLE PRINTOUT PROVIDED BY THE MIAMI 
ILLUSTRATING PERMANENT RECORD OF ICE 
GROWTH HISTORY 
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FIGURE 5 


d(i) 
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Increasing 

Ice 


- variation of MIAMI RESONANCE AS ICE 
ACCRETES ON ITS SURFACE 

(a) Storage oscilloscope trace 

(b) Definition of symbols for use 
with oscilloscope trace 


The microprocessor is programmed to continually monitor 
the change in resonant frequency (indicated by df in Figure 5) 
and convert this change into a dc voltage which is proportional 
to the ice thickness as illustrated in the real time recording 
sho\-m in Figure 2. 

2.2 Performance of the MIAMI 

The MIAMI was tested in the icing research tunnel during 
two weeks in June and one week in August, 1981. Thirty-one tests 
were run in June and twenty-nine in August. 

The purpose of the tests was to determine the ability of 
the MIAMI to measure the thickness of ice accreting on its sur- 
face under conditions simulating those of natural icing. 

The transducer logic assumes that the ice is growing in uni- 
form layers above it. This approximation is true if the trans- 
ducer width w is small relative to major variations in the ice as 
illustrated in Figure 6a. When the ice is non-uniform, as illus- 
trated in figure 6b, the reading produced by the MIAMI represents 
an average value of the ice thickness in the vicinity of the 
transducer. 

For many every day aircraft applications the condition that 
the transducer be small relative to the gross variations in the 
ice layer is satisfied and the MIAMI readings of average thickness 
will closely represent the thickness of the ice layer. For re- 
search applications when the ice layers grow into odd shapes, care 
must be taken that the particular MIAMI chosen for the measurement 
is sufficiently smaller than the gross ice variations in its vicin- 
ity so that the readings are truly representative of the thickness 
of the ice accreting above it. 
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Transducer small relative to variations in ice 



FIGURE 6 - SIZE OF MIAMI TRANSDUCER RELATIVE TO 
GROSS VARIATIONS OF ICE THICKNESS 
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Impact ice under a wide range of tunnel conditions was per- 
mitted to accrete on the MIAMI transducer, which was mounted at 
the stagnation point of a two-dimensional airfoil, as shown in 
Figure 1, and the thickness of the ice was measured by the MIAMI 
and an independent means. The independent means used to measure 
ice thickness during this testing was a calibrated transit located 
in the control room viewing the ice layers accreting on the MIAMI 
surface. The accuracy of this technique, however, is limited to 
relatively thin layers of ice not more than about 125 mils. 

Ice growth on the airfoil under certain cloud conditions was 
non-uniform, being less at the stagnation point than at neighbor- 
ing locations. As the ice grew a time v/as reached when the tran- 
sit could no longer see the ice at the stagnation point, it being 
obscured by the more rapid growth on either side of it. When this 
happened, the transit began to track the more rapidly growing ice, 
or horns, on either side of the stagnation point. This phenomena 
is illustrated in Figure 7 which is one of 18 plots demonstrating 
this phenomena. A more ideal way to test the MIAMI in thick ice 
would be to mount the MIAMI on a rotating cylinder since ice will 
accrete uniformly on this surface. 

Examination of data taken in the icing research tunnel, over- 
whelmingly supports this thesis. The calibration of the MIAMI is 
thus based solely on data gathered during early growth stages when 
the ice is less than about 125 mils (1/8 in) thick and transit 
measurement and MIAMI measurements are in close agreement. 

2.2.1 Galibration of the MIAMI 

Equation 1 represents the theoretically derived relationship 
between ice thickness and shift in resonant frequency as derived 
in Section 3.0. 

I = 1/k In (1-df/df^) (1) 

ci 
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line of transit 
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figure 7 - PEN RECORDING OF MIAMI OUTPUT- ILLUSTRATING HORN GROWTH 


Shift in resonant frequency with ice 


Where df = 
d f* “ 

~ Asymptote of shift in resonant frequency 
with ice. 
k = constant 


I 


1/k ln(l-di/di ) 

cl 


( 2 ) 


Where di 


Shift in frequency index with ice 


^ia = 

k' 


Asymptote of shift in resonant frequency 

with ice 

constant 


It has been found more practical to use Equation 2, however, 
where i, the frequency index or count generated by the computer is 
used in place of frequency since for every frequency index there 
is a corresponding value of frequency f. A measurement of di is 
equivalent to a measurement of df; and di, being continually gene- 
rated by the computer as the ice grows, is readily available from 
computer memory. Curves of measured di vs measured ice thickness 
(e.g.. Figure 8) are readily obtained from the computer printout 
(e.g., Figure 4) and serve as a means for calibrating the MIAMI. 

The calibration process consists of determining values of k' and 
di^ that cause Equation 2 to fit the measured data. Since these 
terms were established from measurements on ice layers up to 125 
mils Equation 2 is, rigorously, applicable only to ice layers less 
than 125 mils. However, since Equation 2 has been verified theo- 
retically, it can be expected to give useful information up to the 
maximum, or dynamic range of the MIAMI transducer. Future tests 
on‘ thicker ice layers are required to verify that the equation is 
applicable to thick layers. 

2.2.2. MIAMI Indications at Temperature N e ar The Freezing 

Point 

In MIAMI tests, run during June and August, a "run" or "test" 
was tacitly considered to be "that period between the time the spray 
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was turned on" and "the time it was turned off." No formal plans 
were made to measure ice growth history after the spray was turned 
off. Accordingly, it was routine procedure to turn the spray on 
and monitor the ice growth history with the MIAMI and a manually 
operated calibrated transit. It was assumed by the test crew and 
the transit operator that when the spray was turned off, the test 
was over and there was no need to continue to measure ice thick- 
ness. Accordingly, the operator of the transit ceased to make 
any further measurements - after the spray was turned off. 

The MIAMI, however, being automatically controlled by the 
micro-processor, did no such thing but continued to monitor and 
record ice thickness and was permitted to run out its preprogrammed 
period of time of approximately 260 seconds. It was during this 
unplanned testing period after the spray was turned off that it 
was noted that when the air temperature was close to freezing, the 
MIAMI indicated an exponential like decay in ice thickness (runs 
1.1 and 2.1). Since there were no independent transit measure- 
ments during this time to compare the MIAMI Indications with, it 
is impossible to establish from the data if there actually was a 
loss in thickness or some other more subtle phenomenon took place 
(at lower temperatures this phenomenon did not manifest itself). 

Various hypotheses explaining this behavior can be made - 
each of which must be verified by additional testing. 

The simplest hypothesis explaining this behavior is that 
there actually was an ice loss at these temperatures. Some evi- 
dence exists (Stallabrass Ref 2) that at temperatures within 
about -4^c of freezing, the adhesive and cohesive forces of ice 
are weaker than they are at lower temperatures, e.g., the ice 


Ref '2 - J.R. Stallabrass and R. D. Price, 

"Ch the Adhesion of Ice to Various Materials" 
National Research Laboratories, July, 1962. 
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structure is not as rigid or is softer than it is at lower tem- 
peratiires and aerodynamic forces acting on the ice structure are 
sufficient to cause it to erode somewhat. As long as the spray 
is on and new ice is forming, the structure can maintain its 
growth. 

A more complicated hypothesis would be the following; When 
the air temperature is at or slightly below freezing, it should 
not be expected that all of the water droplets will freeze instan- 
taneously upon impact. Only a certain fraction will freeze on im- 
pact, forming a rigid ice structure, while the remaining fraction 
of the water droplets will not freeze but will be trapped or 
bound as unfrozen water in the ice structure that has formed. 

There is not sufficient heat transfer during these conditions to 
remove all of the latent heat from the impinging droplets so they 
do not freeze on impact. 

If we confine our attention only to the sample volume within 
the field of the MIAMI, we should expect that within this volume 
some of the bound water in the liquid state is free to migrate or 
drain out of the sample volume, while simultaneously some of the 
bound water will freeze with time, and remain in the sample vol- 
ume. As long as the spray is on new water droplets will be added 
to the sample volume. Some of these will freeze, increasing the 
size of the ice structure within the sample volume and some will 
migrate out of the sample volume. The liquid water migrating out 
of the sample volume being always less than the liquid water being 
added to the sample volume by the spray. The difference between 
the water leaving and the water added is the ice forming within 
the sample volume. 

In summary, what we have within the sample volume in the 
field of the MIAMI, as long as the spray is on, is a continuously 
increasing ice structure containing a certain fraction of un- 
frozen water. If the spray is turned off so that there is no new 
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source of water entering the sample volume some of the liquid 
water will continue to migrate out but some will freeze until a 
point in time is reached when there is no more liquid water in 
the sample volume but only pure ice, or an ice-air mixture. In 
addition, there could be some erosion of the ice structure by 
aerodynamic forces. 

Such an hypothesis will explain the behavior of the MIAMI 
observed when the air temperature was at or near freezing (runs 
1.1 and 1.2). During these experiments, the MIAMI indicated 
linear growth of ice with time, corresponding to the growth in size 
of the ice-water mixture within the sample volume of the MIAMI. 
Since water has a higher dielectric constant than ice, the effec- 
tive dielectric constant of the ice-water mixture within the sample 
volume should be somewhere between the two and higher than that of 
pure ice. When the spray is turned off and the water fraction re- 
duces to zero (by migration of some water out of the sample volume, 
and the freezing of some water that remains in the sample volume), 
the effective dielectric constant will relax or decay, exponen- 
tially, asymptotically approaching the effective dielectric con- 
stant of the ice— air mixture remaining within the sample volume. 

The higher dielectric constant of the ice-water mixture will have 
a greater tuning effect on the MIAMI causing it to indicate thick- 
er ice layers. As the water freezes or is otherwise removed, the 
indications should drop as observed. 

2.3 Measured Dynamic Range of the MIAMI-Model No. 6 

In the 15 runs made with Model No, 6 listed in Table 1, the 
tunnel parameters were set to produce rapid icing so that the ice 
thickness soon exceeded the range of the MIAMI transducer. The 
maximum range of the transducer was taken to be the maximum reli- 
able reading provided by the MIAMI that was not obscured by noise. 
Using this criterion, the average dynamic range recorded was 210 
mils . 
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Model No. 6 operates in the 6 Ghz frequency band. Theory 
(Section 3) indicates that if the operating frequency is lowered 
the dynamic range will increase, (e.g., if the frequency is 
divided by two the dynamic range will double.) Thus, larger 
MIAMI’S should have larger dynamic ranges. Theoretically, dy- 
namic ranges to 1 inch are not unreasonable. 


TABLE NO 1 MIAMI MODEL NO 6, RUNS IN WHICH THE ICE 
THICKNESS EXCEEDED THE DYNAMIC RANGE OF THE MIAMI 



RUN NO 

VEL 

,mph 

1 

LWC. 

*^med 

estimated 
dynamic 
range 
( mils ) 


7.1 

1 50 

0 

1.0 

15 

240 


8.1 

1 50 

-1 0 

1 . 0 

15 

220 


12.1 

200 

0 

1 .0 

15 

1 80 

JUNE 

13.1 

200 

1 0 

1.0 

1 5 

220 

RUNS 

28.1 

1 50 

0 

1 .25 

19 

1 80 


30.1 

1 50 

0 

1 . 8 

19 

220 


32.1 

1 50 

-1 0 

1 .25 

19 

220 


34.1 

1 50 

-1 0 

1 .8 

19 

200 


1 8 

1 50 

20 

1 .8 

19 

220 


20 

200 

20 

1 . 0 

15 

200 


21 

200 

20 

1.5 

20 

220 


27 

1 50 

1 0 

1 .8 

19 

200 


29 

200 

1 0 

1 .0 

15 

1 90 

AUGUST 

30 

200 

10 

1.5 

20 

200 

RUNS 

36 

200 

0 

1 

15 

240 


Average dynamic range = 210 mils 
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2 . 4 Effect of Liquid Water on the MIAMI at Temperatures 

Above Freezing 

As pointed out earlier one of the main features of the MIAMI 
is that it is flush mounted with the surface being monitored. As 
such, it does not interfere with the aerodynamic properties of the 
airfoil. Under these conditions when the surface being monitored 
is placed into an airstream, above freezing, with a high liquid 
water content, the water droplets impinging on the MIAMI will have 
little opportunity to accumulate on its surface. Only small drop- 
lets of water, within the boundary layer of the airstream will re- 
main in the field of the MIAMI. These can be expected to run off 
quickly as their size grows. 

Tests to demonstrate the MIAMI performance under these con- 
ditions were run on June 30, 1981, Run No. 28. The parameters of 
the IRT were set to the following values: 

Temperature 47°F 

Air Velocity 150 mph 

3 

Liquid Water Content 1.8 gm/ cm 

Droplet Size 19 y 

The computer printout and oscillogram indicate only minor, 
negligible, random perturbations of the MIAMI quiescent resonant 
frequency, as might be expected from the above hypothesis. 

It is well known that water is an extremely lossy dielectric 
material and that any significant quantity that is permitted to 
accumulate in the field of the MIAMI will destroy its resonance 
by absorbing the stored energy. This phenomenon was noted in a 
number of unrecorded informal tests performed with water at room 
temperature. One should keep in mind, however, that the phenom- 
enon being discussed here always take place above freezing tem- 
peratures when the MIAMI would not be expected to be operational 
and icing is not expected. 
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3.0 


MIAMI DESIGN THEORY 


T heoretical Basis for the MIAMI 

The theoretical basis for the MIAMI is that the resonant fre 
quency of a resonant surface waveguide is altered or tuned by the 
growth of ice on it. An oscillogram illustrating and verifying 
this phenomena is shown is Figure 5a for MIAMI Model No. 6. 

If the shift in resonance caused by ice of known thickness 
is measured an empirical relation between ice thickness and fre- 
quency shift can be established which, from then on, can be used 
to measure the ice thickness from a measurement of the frequency 
shift. A typical curve of ice thickness vs frequency shift is 
shown in Figure 8 for MIAMI Model No. 6. 

The MIAMI transducer consists of a resonant surface wave- 
guide on which the ice is permitted to accrete and instrumenta- 
tion which continually monitors its resonant frequency and con- 
verts the measured change in frequency to ice thickness by means 
of the empirical equation illustrated in Figure 8 and described 
by Equation 2. The surface waveguide is composed of a layer of 
stable dielectric material such as polyethelene which has approxi- 
mately the same dielectric constant as ice and is placed just 
under the surface on which ice is to be measured as illustrated 
in Figure 9. 

In the ice free condition, the surface waveguide is made 
resonant by making its length equal to an integral number of half 
wavelengths. As layers of ice begin to form on the dielectric 
surface of the waveguide, they have the effect of thickening the 
waveguide and altering its resonant frequency. 

Properties of a Polyethelene Surface Waveguid e 

The properties of a polyethelene surface waveguide may be 
established from a plot of guide wavelength versus frequency for 
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Top View 


FIGURE 9 


i 



- diagram illustrating construction of 

MIAMI TRANSDUCER 
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various thicknesses of the surface waveguide as illustrated, in 
Figure 10, between 1 and 2 Ghz for the TE^mode. This is obtained 
by techniques outlined in Ref. 3. 

Similar plots made over other octave bands such as 2 to 4 
and 4 to 8Ghz, etc., reveal that all plots are identical if the 
scales are related by a normalization factor p so that Figure 10 
may be used to represent any octave band given the normalization 
factor . 


Referring to Figure 10, the surface waveguide of thickness 
d is made resonant by making its length equal to an integral num- 
ber of half wavelengths,. Thus, from Figure 10 when dp = 8 cm. and 
f/p=1.66 Ghz (point a Figure 10), the surface waveguide is reso- 
nant when its length px X = 14.2 cm. This is defined here as the 

g 

quiescent resonant frequency, foo. 


The growth of ice layers, assuming the ice has approximately 
the same dielectric constant as the base dielectric, has the ef- 
fect of increasing the thickness d of the surface waveguide. As 
seen from Figure 10, this has the effect of decreasing the wave- 
length for a given frequency and destroying the resonance. The 
condition of resonance is restored by increasing the wavelength 
by decreasing the frequency until the wavelength is once again 
equal to the guide wavelength that existed in the ice free condi- 
tion. The tuning process is illustrated in Figure 10 as the locus 
of points along the horizontal line of constant . 


Ref. -3 - Feasibility Analysis for a Microwave Deicer 
for Helicopter Rotor Blades USAAMRDL-TR-76-18 , May 1977, 
By: Bertram Magenheim and Frank Hains. 
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First Ord er Approximation of the MIAMI Tuning Curve 

The first order approximation of the tuning curve is obtained 
from this locus described above by making a table of f/p versus 
dp with p X Ag held constant. The approximation is illustrated 
in Figure 11(a). The curve of Figure 11(a) appears to be a de- 
scending exponential and can be very closely approximated by eq 
(3) which is similar to eq (1). (Second approximation Figure 
11(b) is discussed in Section 3.2.) 

X = 1/k In (1-y/y^) (3) 

Where : 

X = d(d p) 

y = d(f/p) 

y = Asymptote of y = 0.5417 

k = constant = 0.268 

3.1 Choosing MIAMI Dimensions 

The dimensions chosen for the experimental MIAMI were based 
on the following considerations; 

• small size 

• widest dynamic range 

• sensitivity 

• available oscillator tuning range 

• oscillator cost 

• ruggedness and reliability. 

The voltage tuned oscillator (VTO) chosen was the Avantek 
Model 8580 that tunes a 12% band from 5.8 Ghz to 6.6 Ghz . This is 
a low cost VTO which, because of its frequency band, would permit 
construction of a small MIAMI transducer. Two designs making use 
of this VTO were prepared, Model No. 5 and 6. Tuning curves for 
these models were plotted on the normalized tuning curve and are 
illustrated in Figure 12. Model no. 5 with p = 3.3 occupies the 

higher sensitivity but lower dynamic range portions of the 


23 



■ Hill II 


■III 



FIGURE 11 - NORMALIZED TUNING CURVE OF MIAMI 
FIRST AND SECOND APPROXIMATIONS 
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1 

lb 

12 
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tuning curve whereas Model no. 6, p = 4 occupies the lower sensi- 
tivity but higher dynamic range portion of the tuning curve. The 
theoretical dimensions of Models 5 and 6 are summarized below. 


Model 
No . p 

Freq . 
Range 

Cut-off 

freq. 

JLength 

a 

Depth 

b 

Width, 

c 

Dyimc 

Range 

Sensi- 

tivity 


Ghz 

Ghz 

inch 

inch 

inch 

inch 

Ghz /in 

5 3.3 

5.8 

5.42 

1.694 

.477 

. 250 

. 263 

5.89 to 


6.6 






3.726 

6 4 

5.92 

3.28 

1.398 

.787 

.25 

.787 

1.96 to 


6.57 .224 

3.2 Construction and Test of Experimental Surface Wave- 
Guides 


Both Model 5 and 6 were constructed and tested. The mea- 
sured resonance curve for Model 5 proved to be non-uniform and un- 
stable possibly because of its proximity to cut-off. Model 6 
proved to be far more tractable; it demonstrated a smooth reso- 
nance curve capable of being tuned by layers of polyethelene tape 
as illustrated in the measurements of Figures 13 to 16. 

The ideal way to tune and calibrate the MIAMI would be to 
place ice layers of precisely known thickness and shape on the 
MIAMI surface; but because of the practical difficulty of doing 
this in the laboratory, at room temperature, the ice layers were 
simulated by adding layers of polyethelene tape of known thickness 
to the MIAMI. Because the tape may be added one layer at a time, 
it could be made to conform to the shape of the MIAMI surface as 
well as serve as a thickness gauge. Table 2 summarizes the tuning 
data taken from Figures 13 to 16. 

The quiescent (no ice) resonant frequency of Model No. 6 
differed from its theoretical value being approximately 6.29 Ghz 
rather than 6.6 Ghz. This is because the end plates of the trans- 
ducer are not perfect short circuits as assumed in the first order 
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FIGURE 13 - PEN RECORDING - MODEL NO . 6 - MEASURED 
SHIFT IN RESONANCE DUE TO LAYERS OF POLYETHELENE TAPE 



-FIGURE 14 - PEN RECORDING - MODEL NO . 6 - MEASURED 
SHIFT IN RESONANCE DUE TO LAYERS OF POLYETHELENE TAPE 
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FIGURE 15 - PEN RECORDING - MODEL NO. 6 - MEASURED SHIFT 
IN RESONANCE DUE TO LAYERS OF POLYETHELENE TAPE 



FIGURE 16 - PEN RECORDING - MODEL NO . 6 - MEASURED SHIFT 
IN RESONANCE DUE TO LAYERS OF POLYETHELENE TAPE 
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Table 2- Tuning Data Ice Simulated by Polyethelene 


(Fig 

No 

■ d(i) 

d(f) 

Ghz 

d{d) j 
mils ■ 

Klein 

BH 

d(f/p) 

Ghz 

d(dp) 

cm 

13 

O 

0 

180 

365 

540 

740 

0 

.042 

.0894 

.142 

.205 

0 

12 

24 i 

48 ! 

96 

0 

.0304 ! 
.061 ‘ 
.122 : 
.244 ■ 

0 

.0105 

.0224 

.0355 

.0513 

0 

.1216 
.244 
. 488 
.975 

u 

0 

0 

0 

0 

0 

0 


21 5 

.0501 

12 

.0304 

.0125 

.1216 

A 

380 

.095 

24 

.061 

.0224 

.244 


469 

.12c 

36 

.091 

.0300 

. 366 


545 

.142 

48 

.122 

.0355 

.488 


680 

.185 

72 

.183 

.046 

.732 


740 

.205 

84 

.21 3 

.0514 

.853 


805 

.227 

96 

.244 

.057 

.975 


880 

.253 

1 08 

.274 

■ 

. 063 

1 . 096 

15 

0 

0 

0 

t ‘ " 

0 

0 

0 


388 

.096 

25 

. 064 

.024 

.254 

□ 

602 

.161 

50 

.127 

.040 

.508 


991 

.282 

1 25 

.31 8 

.071 

1 .27 


— — “ 

.384 

275 

.699 

.096 

2.79 

16 

0 

0 

0 

0 

0 

0 

X 

388 

.096 

25 

.064 

.024 

.254 


600 

.159 

50 

.127 

.0398 

. 508 


625 

.168 

60 

.152 

.0419 

.610 


955 

.279 

— 

1 25 

.318 

.0697 

1 .27 
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theory. They are slightly Inductive and have the effect of elec 
trically lengthening the surface waveguide. Therefore, its qui- 
escent resonant frequency is lower than would be computed from 
consideration of the physical length of the guide alone. The 
adjustment in quiescent resonant frequency is illustrated in Fig 
ure 10 (a to b) and a second approximation of the tuning curve 
obtained due to this adjustment is illustrated in Figure 11(b). 

The portion of the curve occupied by Model No. 6 is indi- 
cated on Figure 11(b). This portion has been replotted in Fig- 

I 

ure 17 fa with coordinates d (f/p) for ordinate and d(dp) for 
absissa. 
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FIGURE 17 


NORMALIZED - THEORETICAL TUNING CURVE 
MODEL NO. 6 - SECOND APPROXIMATION 







3.3 Agreement Between Theory and Experiment Measured 
Tuning; Curves 

Tuning curves for Model No. 6 were measured using simulated 
and real ice as tabulated below. 


Figure 

Symbol 

Dielectric 

Used 

Dielectric 

Constant 

Temp 

Location 

13 

o 

Polyethelene 

2.3 

room 

lab 

14 

A 

II 

2.3 

room 

lab 

15 

CD 

II 

2.3 

room 

lab 

16 

X 

II 

2.3 

room 

lab 

18 


Ice 

3.2 

Freez- 

Wind 





ing 

Tunnel 


The resulting curves (normal 

ized) are illustrated 

in Figure 

18 where they are 

‘ compared with the 

theoretical 

curve (second ap- 

proximation) for 

polyethelene. The 

theoretical 

curve for 

polyethe- 

lene is 

the same 

as Figure 17 fa . 

The tuning 

curve for 

ICE was 


obtained experimentally being the average of 33 test runs in June, 
in the NASA Lewis Icing Research Tunnel (IRT). 


The theoretical curve for ice may be estimated. Since ice 

has a higher dielectric constant than polyethelene (3.2 vs 2.3) it 

will have a greater tuning effect on the MIAMI. An estimate of the 

ice asymptote may be obtained from Figure 10 if the limiting curve 

for ice, A = A //i~"(c = 3»2 for ice) is also plotted and the inter- 
g ^ 

section of and this curve established fa" in Figure 10. The 

ratio of theoretical dynamic ranges being from Figure 10. 


d(f/p)^j 


b-f" 0.47 

a 


= 1.74 


b-f> 0.27 

a 


This establishes d(f/p)al for use in the empirical equation 
eq 3 from which the theoretical curve for ice is established in 
Figure 18 and 17. 
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d(dp) cm 

FIGURE 18 - COMPARISON OF THEORETICAL DATA (SECOND APPROXI 
MATION) AND MEASURED DATA ON NORMALIZED COORDINATES 

MODEL NO. 6 



Where d(f/p)al = Dynamic range for ice 

d(f/p)ap = Dynamic range for polyethelene. 

The difference between the first order approximation for 
polyethelene and measured polyethelene is that the approximation 
assumes the impedance of the end plates of the waveguide does not 
vary as layers of polyethelene are added whereas in reality the 
impedance of the end plates does vary. Similarly, for ice, the 
impedance of the end plates changes as ice is added and since ice 
has a higher dielectric constant than polyethelene the effect is 
even more pronounced. 

3.4 Coupling to the Surface Waveguide (Simplified Theory ) 

The probes used to couple energy into the resonant surface 
waveguide were designed intentionally to provide very light coupling 
so that a high loaded Q with very little radiation would be realized. 
The probes couple energy to the electric fields in the surface wave- 
guide and, therefore, the coupling efficiency is dependent upon the 
relative strength of the electric field at the location of the probe. 
The coupling is zero at frequencies below cut-off of the surface 
waveguide but starts to rise as the operating frequency is increased 
above cut-off. Trying to operate a MIAMI at frequencies only slightly 
above cut-off while yielding very high theoretical sensitivity is im- 
practical because the coupling will be too low. The operating fre- 
quency must be sufficiently above cut-off to provide adequate but 
light coupling. The coupling is also a function of the height of the 
probe relative to the thickness of the dielectric; h/d, as shown in 
Figure 9 (d = doo, when no ice exists h is the height of the probe 
above the bottom of the surface waveguide.) Since the electric field 
at the probe decreases as d increases the coupling will, therefore, 
decrease as ice accretes since this is equivalent to decreasing h/d. 
The probe was placed a quarter wavelength from the end plate since 
this dimension is independent of ice thickness and frequency (guide 
wavelength in a resonant surface waveguide is constant). There is 
no variation in coupling due to this dimension. 
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4.0 


THE MIAMI SYSTEM THEORY 


A block diagram of the MIAMI System is shown in Figure 19 . 
The microprocessor is programmed to generate a repetitive count 
( 1) down from the number 2047 to zero. The count, a number gene- 
rated by the microprocessor, is also called the frequency index. 
The count is converted in the digital to analog converter (DAC) 
to an analog voltage (Vi). This voltage appears as a negative 
going sawtooth the amplitude of which is directly proportional 
to the count. This voltage is transmitted to the MIAMI trans- 
ducer via a long multiconductor cable where it is amplified and 
used to sweep the voltage tuned oscillator (VTO) through its 
frequency range. Each count or number generated by the micro- 
processor corresponds to a different frequency output (f) of the 
VTO. The microwave signal generated by the VTO is then passed 
through the resonant surface waveguide whose resonant frequency 
( f o ) is a function of the ice thickness. The signal is then 
detected by the crystal detector. The voltage out of the crys- 
tal detector (Vo) is then amplified and transmitted via the long 
multiconductor cable to the analog to digital converter where it 
is converted into a digital number (Ai). A peak sensing algo- 
rithm in the microprocessor then determines the frequency index 
(count) that produced a peak output of the resonant surface 
waveguide. This is called the resonant frequency index. 

In the ice free condition, the resonant frequency is called 
the quiescent resonant frequency, foo. The resonant frequency 
index that generated the quiescent resonant frequency is called 
the quiescent frequency index, (ioo). When ice forms on the trans- 
ducer, tuning the surface waveguide, the resonant frequency index 
decreases in value. The difference between the quiescent fre- 
quency index, ioo, and the resonant frequency index, io, is a 
function of ice thickness. 

The functional relationship between ice thickness and fre- 
quency index is established experimentally and stored in the 
microprocessor. The microprocessor then continually monitors the 
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Ice growth 
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FIGURE NO 19 - MIAMI SYSTEM BLOCK DIAGRAM 
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change in frequency index and using the stored function calcu- 
lates the ice thickness. The microprocessor then computes the ac 
cretion rate of the ice by calculating the first time derivative 
of ice thickness. These values are then available for subsequent 
display . 

4^1 Th eore tic al ly Derived Relation Between Ice Thickness 

and Frequency Index 

The relationship between frequency index and ice thickness 
may be obtained functionally from Figure No. 19 provided the 
transfer constants of the surface waveguide as a function of ice, 
the VTO, the VTO driver and DAC are known. All of these trans- 
fer characteristics can be measured and derivation of their 
values and functional relationships appear in Appendix No.l 

From eq (1) the ice thickness is a function of the shift 
in resonant frequency. 

I = 1/k ln(l-df/df ) (4) 

The VTO tuning curve is represented by 

f = f - f e®^ (5) 

am 

and the transfer characteristics of the DAC and VTO driver are 

V = mi + b (6) 

and 

dv = m di ( 7 ) 

from which it can be shown that df is related to di by 

df = (8) 
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Thus, the ice as a function of di can be obtained from eq 4 and 
8 as 


I = 1/k In 


_ gv [ -gm di T ^ 
T f e^ oo (e ^ -1) 
1— m 

df 


(9) 


Parameters in eqs 4 to 9 are defined in the following table: 


Where 


Poly. 

Theoretical 

Poly. 

Measured 

Ice 

Theoretical 

Ice 

Measured 

df 

frequency shift 





di 

index shift 





dfa = 

frequiency asymptote 

CM 

00 

O 

0.354 

1.4268 

.573 

k 

constant 

~2 . 5 

-11.45 

-2.5 

-9 . 51 

foo 

quiescent frequency 


6.29 Ghz 



fm = 

VTO parameter 


1 . 30 



g 

VTO parameter 


-0.0735 



m 

slope of VTO driver 


0.0058 



Voo = 

quiescent VTO voltage 


11 . 68 



e 

base of natural log 






Figure 20 was prepared to illustrate the agreement between the 

theoretically derived eq 9 and measurements. Four plots of eq 9 

are shown for ; 

Ice measured 
Ice theoretical 
Polyethelene measured 
Polyethelene theoretical. 

Parameters of the MIAMI components tabulated above do not 
change in any significant manner compared to the shift in resonant 
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FIGURE 20 


COMPARISON OF THEORETICAL AND EXPERIMENTAL 
DATA OBTAINED WITH MATHEMATICAL MODEL OF 
MIAMI SYSTEM 


39 



frequency index. The microprocessor monitors changes in index (di) 
not absolute values of index (i) so that the effect of changes in 
foo or ioo with ambient temperature are minimized. All component 
parameters tabulated above were established experimentally. _ 

4.2 Variation of Dynamic Range with MIAMI Di m ensions 

The dynamic range of the MIAMI system is limited primarily 
by the design of the MIAMI transducer. The dynamic range for a 
particular transducer can be established from its tuning curve. 

The tuning curve for Model No. 6, as measured during tests performed 
in June 1981, is illustrated in Figure 18a on normalized coordinates. 
The practical normalized dynamic range is indicated as ”rp" on Fig- 
ure 18a (r is the dynamic range, p is the normalization factor). 

The dynamic range is defined here as the ice thickness indicated by 
the MIAMI when the resonant frequency of the MIAMI has shifted to 
a value equal to 80% of the frequency asymptote (in Figure 18a 
the value of normalized dynamic range, (rp) is 2.0 cm when d(f/p) 
is 0.180 Ghz which is approximately 80% of the asymptote of d(f/p) 
which is about 0.225). 

The dynamic range thus defined is useful in that it is a num- 
ber representing the order of magnitude of ice thickness that can 
be accurately read by the MIAMI. The arbitrarily chosen 80% cri- 
terion may be more conservative than necessary since it may be 
possible to obtain accurate readings from the MIAMI when the fre- 
quency shift is closer to the asymptote. 

Thus, using a conservative normalized dynamic range rp of 
2.0 cm, the following table shows how the absolute value of dynam- 
ic range, r, varies as the normalization factor, p, is varied. 
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Variation in Dynamic Range, r, with Normalization 
Factor, p, established from MIAMI Model No. 6. 


Normalization 

Factor 

Dynamic 

Range 

if rp =2cm 

Quiescent 

Resonant 

Frequency 

Thickness of 
Surface 
Waveguide 


P 

r 

foo 

doo 




cm_._ in._ 

Ghz 

_ in . .. 

cm 


0.5 

4 1.57 

0.787 

6.29 

16 


1 

2 . 787 

1.573 

3.149 

8 


2 

1 .393 

3.145 

1.574 

4 


4 

0.5 .196 

6 . 29 

0.787 

2 

'Measured 


Physically, the dynamic range is limited by the following con- 
siderations. The low end of the tuning range becomes useful only 
when the operating frequency is far enough above cut-off for the 
coupling efficiency to be sufficiently high to couple useful 
energy into the surface waveguide. The high end of the tuning 
curve becomes useless when the asymptote is approached and a 
small change in ice thickness produces only an unmeasurable change 
in resonant frequency. 

This table shows that the dynamic range of the MIAMI may be 
adjusted upwards or downwards to accomodate the particular ice mea- 
surement at hand. 

This table shows that when the normalization factor, p, is 
halved, the dynamic range doubles. The normalization factor used 
in the design of Model No. 6 was p=4, and ice thickness readings 
to at least 0.200 (0.196) were recorded. The above table shows 
that if the normalization factor is dropped from 4 to 1, the dy- 
namic range will increase from 0.196 inches to 0.787 inches with a 
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corresponding increase in MIAMI depth doo from 0.787 inches to 
3.149 inches. The quiescent frequency foo would aecrease from 
6.29 Ghz to 1.573 Ghz . 
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5.0 RECOMMENDATION 


It is recommended that a program be initiated that will 
improve, miniaturize and ruggedize the Microwave Ice Accretion 
Measurement Instrument (MIAMI), so that it can be mounted on a 
rotating cylinder and, in a later program, on fixed or rotary 
wing aircraft for research purposes. It will answer questions 
raised concerning the dynamic range, and the effects of unfrozen 
water and other liquids on MIAMI performance. 

Calibration and the Effects of Unfrozen Water 

An ideal way to calibrate the MIAMI is to mount it on a 
rotating cylinder since ice will accrete uniformly in a predict- 
able and easily measured manner on this surface. 

The installation will be used for the following purposes; 

1. Calibration of the Miami. 

2. Determination of the MIAMI dynamic range. 

3. Establish the relationship between dynamic 
range, operating frequency, and MIAMI size. 

4. Determine the effects of unfrozen water on 
MIAMI ' s performance . 

5. Measure the unfrozen water content of the 
ice . 

6. Measure the variation in unfrozen water 
content with time (a transient phenomena). 

7. Improve system software. 

Mounting of the MIAMI on Aircraft Surfaces 

Once the improved MIAMI transducer is calibrated it may, at 
some later time, be transferred to either a fixed or rotary wing 
aircraft for the study of the effects of naturally occurring ice 
on aircraft performance and general research into the icing phe- 
nomena and cloud research. 
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Many dangerous icing encounters are transient in nature and 
disappear before they can be measured. With an installed MIAMI, 
a complete history of every icing encounter can be objectively 
recorded . 
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6.0 


LIST OF SYMBOLS 


Symbol 

di 

w 

k, k' 

f 

df 

I 

dl/dt 

i or j 

di or dj 

di or dia 
a 

foo 
f o 

ioo or joo 
^o or jo 
T 

d , 
med 

LWC 

VEL 

d 

doo 

^go 

P 

TE 

Ghz 

g 


Definition 

Shift in index 

Width of MIAMI transducer 

A constant used in empirical 
calibration curve 

Frequency 

Change in resonant frequency 

Asymptote of change in resonant 
frequency 

Ice thickness 

Accretion rate 

Frequency index (used interchangeably) 

Change in frequency index 

Asymptote of change in resonant 
frequency 

Quiescent resonant frequency 
Resonant frequency 
Quiescent frequency index 
Resonant frequency index 
Temperature 
Droplet size 

Liquid water content 
Airspeed 

Thickness of surface waveguide 

Thickness of surface waveguide used in 
transducer which is resonant at foo 

Free space wavelength 

Guide wavelength of surface waveguide 

Guide wavelength of quiescent resonant 
frequency 

Normalization factor 
Transverse electric 
Giga - Hertz 
Parameter of VTO 
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Symbol 

m 

Voo 

e 

fm 

r 

0 - 


Def inition 

Parameter of VTO 

Quiescent VTO voltage 

Base of natural logarithm 

Parameter of VTO 

Dynamic range of MIAMI 

Quality factor of a resonant device 
(ratio of energy stored to energy 
dissipated) 


46 



APPENDIX 1 


1.0 TRANSFER CHARACTERISTICS OF PAG AND THE VTO DRIVER 

The output voltage of the VTO driver amplifier is a linear 
function of the frequency index J (j is used interchangeably with i.) 
V = mj + b (9) 

where m and b have been established experimentally for Model No. 6 

m = 0.0058 volts per index 

b == 2.75 volts. 
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2.0 


TRANSFER CHARACTERISTICS OF VTO 8580 MEASURED 


f = 


f 


a 



( 10 ) 


Where 

f = frequency output of VTO in GHz 
f = asymptote of VTO tuning curve = 6.85 GHz 

3 . 

f^ = y intercept of VTO tuning curve 

when plotted on semi log paper = 1.3 

g = slope of VTO tuning curve 

when plotted on semi log paper = -0.0735 


note : 

Power output of VTO 8580 
= 5mw + 1.5 db 

When f = 6.6 Ghz ^ V = 24.0 volts 
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3.0 FREQUENCY TRANSFER CHARACTERISTICS OF VTO (DERIVATION) 

3.0,1 The following diagram represents the theoretical timing 
curve of a varactor tuned oscillator 



when Y = CO 

e^^= 0 (since g is negative) 
f = f 

a 

when V = 0 



solving eq(11) for g 

In (f - f)/f 
a m 

g = 

V 


( 12 ) 


rewriting eq (11) 

(f - f) = f e^^ 
a m 


(13) 


In (f - f 
a 

equation (14) is 
paper with slope 
in the following 


) = gv + In f 
^ m 


(14) 


the equation of a straight line on semi log 
= g and y intercept = In f^ as illustrated 
diagram . 
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To fit the experimental data to the equation, eq (ll). 

1 . Estimate the value of the asymptote f by plotting experi- 

3 . 

mental tuning data on rectangular coordinate paper as in Fig. 21. 


2. With initial value of asymptote f plot eq(14)'on semi log 

sl 

paper Figure 22, The plot should be a straight line. If all 
points fall on the straight line, the initial value of f chosen 

d. 

was a good approximation. If values do not fall on the line, 
choose values for f by trial and error until all points fall on 
the line as illustrated in Figure 22. The highest frequency point 
closest to the asymptote will be the most sensitive to varia- 
tions in f when tli.s point falls on the line the approximation is 
very close. The following is a table summarizing the successive 
estimates of fa. 
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3 . 


The final constants of the previous example are 
f = 6.85 Ghz, f = 1.5 Ghz, g = -0.0735 

These constants when used in equation 11 yield the curve of 

Figure 21, illustrating an excellent fit to the experimental 

data. (Later measurements show that f = 1.3 provides a bet- 

m 

ter fit to the VTO used in Model No. 6.) (Note: fm = fa— fb, 
fa was just estimated and fb is read directly from Figure 21.) 
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FIGURE 21 - VARACTOR TUNED OSCILLATOR TUNING CURVE 
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4.0 MEASURED CHANGE IN FREQUENCY (df) WITH CHANGE IN FREQUENCY 

INDEX (di) 

The change in frequency due to the applied frequency index 
was measured in the following manner 

Known thicknesses of polyethelene were applied to the MIAMI 
transducer and the frequency shift produced was measured in the 
laboratory at room temperature. Data is presented in Figures IS- 
IS and summarized in Table 2 and is plotted in Figure 23. Knov/n 
thicknesses of polyethelene tape were also applied to MIAMI trans- 
ducer and the change in applied frequency index measured, and are 
plotted in Figure 24. 

From Figure 23 and 24, the change in frequency produced by 
a change in index is established and is graphed in Figure 25. 

The change in frequency due to change in index may also be 
established from the VTO tuning curve given the relationship be- 
tween drive voltage and frequency index and the constants of the 
VTO as presented in Equation 10, Good agreement between eq 10 
and Figure 25 is achieved when fm = 1.3, g = 0.035m = 0.0058 and 
voo = 11.68 volts. 
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I (nils) 

FIGURE 23 - RESONANT FREQUENCY SHIFT (MEASURED) DUE 

TO LAYERS OF POLYETHELENE ON MIAMI TRANSDUCER 
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FIGURE 24 - CHANGE IN FREQUENCY INDEX (MEASURED) 

DUE TO POLYETHELENE ON MIAMI TRANSDUCER 
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Ghz 



FIGURE 25 - FREQUENCY AND CHANGE IN F’REQUENCY AS A 
FUNCTION OF CHANGE IN FREQUENCY INDEX - 
MIAMI MODEL NO. 6 FROM JUNE TESTS 
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5,0 EQUIVALENT CIRCUIT THEORY OF MIAMI TRANSDUCER 

A physical circuit of the MIAMI transducer appears in Fig- 
ure 26 with a layer of ice that is smaller than the base dielectric 
doo . The thickness of the base dielectric doo is designed so that 
in the absence of ice the surface waveguide between terminals 1 
and 2 can support the TEl mode surface wave and therefore can be 
considered a "uniform transmission line" of thickness doo. When 
an ice layer of thickness I, with I always less than doo, is 
added the modified waveguide between terminals 1 and 2 can still 
support the TEl surface wave, however, since the ice layer of 
thickness I is always less than doo, the ice outside terminals 1 
and 2 cannot support the TEl mode. This ice appears as a wave- 
guide beyond cut-off to the surface waveguide between terminals 
1 and 2, and presents an inductive load to the waveguide at termi- 
nals 1 and 2 as shown in Figure 26b. Even in the absence of ice, 
the end plates at terminals 1 and 2 are not perfect short circuits 
but also present an inductive load to the waveguide. Thus, the 
inductance shown in Figure 26b has contributions from the end 
plates as well as the ice layers. 

From transmission line considerations (Smith Chart) a short 
length of transmission line is inductive so that the inductance 
of Figure 26b may be replaced by two short transmission lines of 
length 1*. The waveguide is resonant when the equivalent length 
of transmission line is equal to an integral number of half guide 
wavelengths . 

The theoretical tuning curve for a polyethelene surface 
waveguide was derived assuming the inductances presented to the 
transmission line were fixed. Since these vary also with ice 
thickness, we should expect the measured tuning curve to deviate 
somewhat from the theoretical curve. 
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FIGURE 26 ~ EQUIVALENT CIRCUIT - MIAMI TRANSDUCER 



APPENDIX 2 


SYSTEM DOCUMENTATION 


1.0 THE MICROPROCESSOR 


The general structure of the MIAMI, including the micropro- 
cessor is shown in Figure 27. In operation, the computer generates 
a number J which is presented to the DAC (a)l at its memory loca- 
tion 08ff0H. The DAC converts the number to a current that is 
linearly related to J by I = J + which ranges between 4 and 

20 ma for J = 1 to 2047. This current is then transmitted to the 

transducer via a long cable where it is impressed across the input 
resistor of the VTO driver where it is converted into a voltage and 
is amplified by the VTO to the 0 to 30 volt range for driving the 
VTO. The VTO in response to the tuning voltage oscillates in the 
5.8 to 6.6 Ghz band, its exact frequency being a function of the 
tuning voltage. 

The signal is then passed through the resonant surface wave- 
guide and then detected by the crystal detector resulting in a 
negative voltage between 0 and -0.5 volts. This voltage is then 
amplified by the detector amplifier and transmitted via the long 
cable to the ADC boards in the computer. The returned voltage Aj 
is presented simultaneously to two ADC converters, one with unity 
gain and the other with a gain of 3 for use when the signal is 
weak. After conversion by the ADC’s, one or another of the two 
resulting numbers is read by the computer and represents the re- 
sponse of the resonant surface waveguide (plus its ice load) to 
the applied VTO frequency f (J). 

The computer then emits another value of J and processes 
the returns, looking for a maximum Aj . When a maximum Aj is 
found, the value of J, which produced it, ML, is the "resonant 
frequency index." The quiescent frequency f^ (no ice) produces 
a frequency index CA or calibration index.* A function if CA-ML 
= DX is then related to the ice thickness at the transducer. 


*As of 9/9/81, if CA=1505, foo=6.29 Ghz, MX=1465 
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When ML is found, it is used by the computer to determine 
the range of values of J to be presented next, J1 to J2, and the 
process of searching for the peak between J1 and J2 is repeated, J 
again ranging downward in unit steps from J1 to J2. The sampling 
rate is 6000 frequencies/sec., every one half second for the first 
30 seconds, and every second thereafter. At the completion of a 
sweep, the latest values of DX, the time of the on-board clock T, 
and ice thickness 0, and the amplitude of the peak (number pro- 
portional to its voltage) MX, are stored in an array in computer 
memory. At the same time, one of three possible pairs of values 
is presented to the oscilloscope and plotter input terminals, and 
ice thickness and accretion rate are converted to decimal (from 
the computer's binary) and presented to the two panel digital dis- 
plays ( PDD ' s ) . 


Data available for online oscilloscope/plotter output are 
the pairs: (time, 0), (time, ML), (J, A j ) , where the first term ap- 

pears as the x-coordinate , and second as y. While the oscillo- 
scope and plotter are in general operated in parallel, it is pos- 
sible to read J, A^ , on the oscilloscope while plotting any of 
the other pairs at the same time. 

The program as stored on EPROM does not provide for alternate 
displays, nor different kinds of plots, outputs are J, A , during 

yj 

test, MX, ML after calibration and 0, T after icing starts. For 
alternative outputs, the program must be read in from cassette 
tape and proper modifications made in the code. Alternate post- 
run plots are also coded (but not linked) in the taped program. 

Interrupts (to switch functions of the computer) are pre- 
sented via key strokes on the printer during execution. For cali- 
bration of ice thickness as seen by the MIAMI with ice thickness 
as seen by a transit, a "Mark!" stroke, entered via the keyboard 
as a CR (carriage return), records the time at which the event 
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occurred. (See operating instructions.) At the completion of a 
run, the data stored are printed, and other plots can be made, in 
pairs, from the set T (time), 9, 9, ML, MX, where 9 is accretion 
rate (d9/dt). Then the data are erased and the program reset for 
the next run. 

The two most frequently used subroutines, "PEAK'* (which sweeps 
the range of frequencies and determines the peak) and "HEXDEC" 
(which converts binary to decimal for PDD display) are coded in 
assembly language. The rest of the program is BASIC. 
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1.1 


Peak Following Algorithm 



J is decremented from J1 to J2. A is amplitude of 
return signal . 

MX is highest value of A found. 

MJ is value of J at which MX is first found. 

MK is value of J at which MX is last found. (May be equal to MJ ) 
ML = (MJ + MK)/2. 

CA is ice-free (quiescent resonant frequency) index. 

Initial value of J1 = CA + 32, J2 = CA -32. 

New values of J1 and J2, after sweep is completed, 

are J1 = ML-48 , J2 = ML+48, unless 

MK-J2 < 3 in which case J1 = J2 + 16, J2 = Jl-96 

or 

J1 - ML < 3 in which case J2 = J1 - 16, J1 = J2 +96 
or 

ML > CA in which case J1 = CA +8 , J2 = CA - 88 
or 

MX < 500 or OS > 8 in which case J1 = CA, J2 = 1 

where OS is amplification flag. If OS = 8, 9,10 or 1 1 
then A is amplified by 3 and sweep covers entire range. 
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If MX < 500 and OS < 8 then OS = OS + 8 
If MX > 2000 and OS > 7 then OS - OS - 8 

This provides 'hysteresis' to prevent oscillation. During 
test (before calibration) J1 = 2000, J2 = 1 


1.2 Averaging Subroutine 


During calibration , to the nearest integer 

30 

CA= 1/30 ^ MLj 

1 = 1 

The ice trigger is 


IT = CA - 2*MIN(MLj) 

When waiting for ice if in any sweep 
CA-ML> IT 


the assumption is that icing has begun. 

Ice thickness smoothed TT and accretion rate TK are determined 
from the six most recent readings of TH=f(CA-ML) and fitted by 
a least squares straight line. 
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2.0 THE COMPUTER PROGRAM 


2 . 1 Flow Charts 


HIAHI CoHputer ProgrM 


INVKATION 

INITIALIZE 


TEST 


calibrate 


SET UP 
DISPLAY 


WAIT FOR 
ICE 



FIGURE 28 - MIAMI COMPUTER FLOW CHART 
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COMPUTER FLOW CHART 



FIGURE 30-MIAMI COMPUTER FLOW CHART 
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COMPUTER FLOW CHART 


ASSEMBLY PROGRAM "PEAK" 


set up working storage 
find and store return address 
clear status words, both boards 
select input channel, both boards 


A4 


A2 

A1 

A3 

A5 

A6 

A7 


move J1 to J 

load -999 into A (R5) 

load -1000 into MX (r6) 

move J — ♦ DAC A output (next cycle) 

start conversion both boards this cycle 

compare last value of A with MX 

if A < MX* go to A1 

if A ==mx go to A2 

move A to MX 

move J to MJ 

move J to MK 

compare OS flag with 8 bit (gain) 
if low gain go to A5 

if conversion is not complete, go to A3 * 
move input data ('A') to high gain 
jump to A7 

if conversion is not complete, go to A5 * 

move input data ('A') to low gain board 

if 0S<’>1 (display A,J) go to A7 

decrement J output A,J to ADC 

if J < J2 go to A4 

calculate (MJ+MK)/2 = ML 

return 


( spin here until conversion is complete) 


FIGURE 31 
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2.2 Main Control Program, Definitions 


T(3) 

V(241,4) 
U(30, 4) 

C2 ,C3,C4, 
C5, 


storage for time 
data table for time, MX, TH, DX 
transit reading tabletime, DX, TRANSIT ICE, 
MIAMI ICE 

coefficients of fifth order polynomial 

0= CCj^DX^ 

i = 1 


$B 

CA 

ML 

DX 

CU 

U 

V 

W 

MX 

MK 

MJ 

MWD(W+6 ) 
TT 

DEF(DX) 

U1 

U2 

VK 

S 

J1 

J2 

OS 


Z1 , Z2 
A 
SM 
MM 


token for beeper 

calibration index nominally 1770 

resonant frequency index 

CA-ML 

row in v table at which to start smoothing 
memory address (hex) of ADC(l) data block 
memory address (hex) of ADC ( 2 ) data block 
memory address (hex) of workspace for assembly 
'PEAK' 

maximum value of Aj 

last (lowest) value of J for which A=MX 
first (highest) value of j for which A=MX 
increment between steps in "PEAK" subroutine 
smoothed ice thickness 
DXd(C2+DX(C3+DX(C4+DX(C5+DX*C6) ) ) ) 
row index of U table 
flag, U table overflow 
row index of V table 
step in J same as MWD(W+6) 
initial (start) value of sweep frequency 
index 

final (end) value of sweep frequency index 
amplification and display flag 

8.0 no display 

9.1 amplitude j 

AH peak amplitude J 
CH ice thickness, time 
if OS > 7 amplification is 3 
if os < 8 amplification is 1 
truth table 

OS 0-500 500-2000 2000 

■'8 OS — 8 

>7 - - 0S-»<8 

arguments for 'HEXDEC" 

keyboard interrupt character 

sum of readings of MK during calibration 

minimum reading of MK during calibration 
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2.2 Definitions (cont*d) 


CT 

IT 

TN 

TM 

T 

TD 

TH 

TM 

linel075 


$A 

M 


TT 

TK 

S1,S2,S3, 

S4,S5 

DN 

Ml 

A1 


J 

K1 ,K2 

f(J) 

ML 


CA 

DX 

J1 

J2 

MX 

T,t 

CR 

0, TH 
TT 


if CU=13, storage interval is 13/25 seconds 
if CU=25, storage interval is 25/25 seconds 
ice trigger 

time since clock started, equals moment 
ice starts to form 

equates TN, first storage interval in 25 
of a second (1 tick = 40 ms.) 
time since start of icing 
time since last storage 
ice thickness 
time now defined in 1050 
ASCII CHAR 

13 carriage return 

65 A 

66 B 

67 C 

68 D 

69 E 

input string, printed, not saved 
transit coefficient. Since this number 
(.93) is stored on EPROM if a different 
value is used U table values of observed 
ice must be adjusted manually- 
smoothed value of ice thickness 
smoothed accretion 

terms in least squares fit 
denominator, if zero start at next point 
accretion rate mils/second 
intercept 

SYSTEM DESIGN NOTATION 

index number 0<J<2048 
constants: local to expression 

amplitude of return signal befo.re ADC 
frequency generated by VTO 
ML=(MJ+MK) /2 

quiescent resonant frequency 

quiescent frequency index (calibration value) 
CA-ML 

first value of J presented during sweep 
last value of J presented during sweep 
J1>J2 

amplitude of peak (maximum value of A found) 

time in seconds 

carriage return 

ice thickness 

ice thickness smoothed 
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SAVE 2 
NEU 

TM990 BASIC REV B.1.10 
♦READY 
LOAD 05 00 OH 
LI S 

10 REM ♦♦♦♦MIAMI CONTROL PROGRAM#^^^ 

11 REM IDEAL RESEARCH# INC. # ROCKVILLE# MB.’ 

12 REM AUTHOR r JAMES K.- ROCKS 

13 REM DEVELOPED UNDER NASA CONTRACT# 1981 

14 REM ALL RIGHTS RESERVED 

20 DIM Tt33#Vt241#4J#U£30#41 

25 TIME 0#0#0 

30 TIME ITfOI 

40 C2=0.- 0426936341 C3»-*0.‘ 00009563097: t C4»3.'50247865E-‘07 

42 C5=-3.-49923332E-10:: C6=l .-25303637E-13: : SB=X7JsO:: DX=0:: CU»7 

44 CA=1770:: U=03FF0H:: V»09FF0H:: W-OC700H:: MK=0:: MJ*0:: ML»0 

46 MWDW+63*-!:: TT=0 

50 DEF FNDr8J*£^<e2+t^<C3+t»«;C4>r^<C5+f^C6>»> 

405 DATA 0C700M#0e003M#0420H#eC000H 

410 DATA 020eH#016H#0A00DH#0C390M#04E0H#03FF4H 

420 DATA O4E0H# 09FF4M# 0200H# OlH# 0C300H# 03FF8H# 0C300H# 09FF8H# OCOOIH# 0205M 

430 DATA 0FC19M# 0206H# 0FC18H# 0C800H# 08FF0H# 0720M# 08FFAM# O720M# 09FFAH# 03185H 

44 0 DATA O11O4H#O1302H#OC1 35H # OC 1 C OH# 0C2 0 OH# 02 1 2 OH # OC 034M # 0 1 6 06M # 056 OH 

450 DATA 09FFCM# 01 IFDH# OCieOM# 09FFEH# 01 005H# 0560H# 08FFCM# 01 IFDM# 0C160H# 08FFEM 

460 DATA 02120M#OC07EH#01604H# 0C3O5H# 09FF2M# 0C800H# 09FF0H# OA003H# 03030H 

470 DATA 015DCH#OC247H#OA243M#0819M#0330M#01H#02H#04M#08H 

500 FOR l=0ee00H TO 0C084M STEP 2:: READ A:: MWDE1J=A:: NEXT 1 

510 DATA 0C720H# OCl 03H# 0420H# OCl OOH# 0200H# 016H# OAOODH# 0C390H# OCOAOM 

512 DATA Oe7FAM#O4C0H#O2O4H#O4H#O2O5M#OAH#OB4OH#OeO32H#O13OCH 

514 DATA OeOC2M# 04C2H# 03C85M# 0E003H# 0604H# 016F7H# 0B40H# 0C042H# 01 OOOM 

516 BATA 01000M#01000M#e330M#OB40H#0604M#016FDH#010F7M 

520 FOR l=0ei00M TO 0C142H STEP 2:: READ A:: NWDfIJ*A:: NEXT 1 

600 FOR 1=1 TO 240:: FOR J=1 TO 4:: V£I#JJ=0:: NEXT J:: NEXT I 

610 FOR 1=1 TO 30:: FDR J=1 TO 4:: U£I#JJ=0:: NEXT J:: NEXT 1 

620 Ul=i:: U2=0:: VK=1 

700 S=-i:: Jl=2000:: J2=l:: DS=1 

710 GOSUB 5010 

730 Z1=MX:: 22=0220H:: GOSUB 4000:: 21=NL:: 22=02A0H:: GOSUB 4000 

750 A=NKYE0I:: IF A 067 THEN GOTO 710 
310 03=1:: SM=0:: MM=2000:: CT=13 

320 Jl=ML+32:: J2=ML-32 

330 FDR 1=1 TO 30:: GOSUB 5010:: SM=SM+ML 
340 IF MK<MM THEN MM=MK 

350 NEXT 1:: CA=INP£SM-'30+0.-51:: IT=2^<CA-MM> 

360 VE1#1I=0:: V£1#23=MX:: VE1#33=CA:: V£1#4J=0:: VK=2 

905 21=CA:: 22=0220H:: GOSUB 4000:: 21=IT:: Z2=02A0H:: GOSUB 4000 

910 MUD£09FF0HJ=20O0:: MWD E09FF2H3-2000 

915 0S=4:: MUD £09FF0Hl=O: : MWD E09FF2H3=0 

920 Jl=CA+64:: J2=CA-64 

930 GOSUB 5010 

935 A=NKY£0J:: IF A-67 THEN SOTO 810 
940 IF CA^ML<1T THEN GOTO 910 
960 DS=4:: TN=TICE0J:: PRINT fB:: TM=TN 

970 J2=ML-48:: J1=ML^48 

930 GOSUB 5010 

990 IF OS <3 AND MX>500 THEN GOTO 1010 

995 IF 0S>7 AND MX <2000 THEN GOTO 1007 

1000 IF 0S<8 AND MX<500 THEN DS=OS-»«:: GOTO 1007 

1005 IF DS>7 AND MX>2000 THEN 03=03-8:: GOTO 1010 

1007 J1*CA:: J2=l:: GOTO 1030 

1010 IF ABS£MK-J23<3 THEN J1=J2+16?: J2»J2-30:: GOTO 1030 
1015 IF ABSEJl-MJIO THEN J2*J1-16:: J1»J1^80:: GOTO 1030 
1017 Jl»ML+48:: J2*ML-48 

1030 T=T1C£TNT:: TD*T1C£TM3 
1040 IF TD<CT THEN GOTO 980 
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1845 IF BX<8 THEM TM-8?? Jl=Cft+0fr J2*CH-88 

1 846 TM*FND CDX J 

2858 TM=TICf8Jr? VEVK# ll»T/25tr V tVKr2J»MXs ? V EVKir3J»Cft-MLs f VfVRir4J=TH 
1875 ft*N(CYE8J.-? IF H013 THEM SOTO 1882 

1838 PRINT -.-irr U tUl if 1 J=T#'8.- 84r t U1*U1+1?? IF Ul>38 THEM Ul»38?? U2=02+l 

1832 IF fl=66 THEN SOTB 1588 

1834 IF VK>CU THEN SDSUB 7888 

1836 IF VK»CU THEN SQSUB 7288 

1898 VK=VK-M 

1188 IF VK<68 THEN OTTn 988 
1118 CT=25 

1128 IF VK<241 THEN 6DTD 938 

1580 INPUT -RUM NUMBER? INPUT -TUNNEL PRHNS? VEL*->fH? 

1528 INPUT - TEMPi[F>=*-?»ft» ?? INPUT - BRDPLET SI2E=-»lft» 

1530 INPUT - WHTER CONTENT*- If H 
1540 INPUT -COMMENTS? -»fP 
1545 N=0.-93 

1558 PRINT -TRHNSIT COEFFICIENT? -IN?? PRINT ?? PRINT 
1340 J=1 

1350 FOR 1=1 TO VK 

I860 IF UtJy lJ>VtI)f II THEN SOTO 1890 

1370 UEJ5r4I=VtIi-4I?? UEJjr2J=VEI)r3J?? UEJ>3I*M^J^28?? J=J-M 
1330 IF J>=U1 THEN 6DT0 1988 
1390 NEXT I 

1900 PRINT - TIME D-INDEX TRHNSIT ICE NIHMI ICE- 

1918 FDR 1 = 1 TO Ul?? FDR J=1 TO 4?? PRINT ??-SSSS,99 -UElifJlJ 

1920 NEXT J?? PRINT ?? NEXT I 

1930 PRINT -TlME-^-RESDNftNT-!f-FREOUENCY-!f-TMICKNESS-f-ftCCRETIDN- 
1941 PRINT -SECS-jf-ftMPLITUDE-jf- SHIFT- «r -MILS OF ICE- .f -RftTEj»MILS<^NIN- 
1945 VH=Vk 

1950 FOR 1=1 TO VH?? PRINT V El j» 1 J f V El »2J sf V El jf 31 jr 

1955 PRINT c-SSSSS.-9-VEI)f4II 

1956 T=VEIifl3 

I960 IF KCU THEN PRINT ?? GOTO 1966 

1962 VK=I?? IF I>CU THEN 60SUB 7880?? GOTO 1965 

1964 IF I=eU THEN GOSUB 7280 

1965 PRINT - '-lJ?-‘SSS.-99 -TT»- -?«-SSS.-99-Tk 

1966 H=NKYE8I?? IF H=66 THEN GOTO 1972 
1970 NEXT I 

1972 PRINT ?? PRINT 
1975 VK=VH 

1930 PRINT -ENTER R FDR NEXT TESTED TO REDO REPORT# E TO REBO REPORT WITHOUT H 
DR?- 

1931 PRINT ?? PRINT ?? PRINT ?? PRINT ?? PRINT 
1982 R=NKYE83?? IF 3=8 THEN GOTO 1982 

1935 IF H<65 DR ft>69 THEN GOTO 1938 

1990 ON H-64 THEN GOTO 688# 1930# 1938# 1508# 1348 

4 000 MWBE0C7FftHJ=Zl?? CHLL '‘HEXBEC" # 8C1 84H? ? BftSE 22?? CRF E8I=MWB E8C728HI 
4010 RETURN 

5010 MWDEW+2J=J1?? MWBfW+4I=J2?? NWB f W+8 J =03 
5020 CHLL -PEftk"#eC804M 

5030 MX=MWBEW>12J?? ML=MWB EW+18J ? ? MJ=MWB EW+14J ? ? NK=MWBEW+16J 
5040 RETURN 

7000 S1 = S1+VEVK#1I^VEVK#4J-VEVK.-6#1I»VEVK-6#4I 

7010 S2=S2+VEVK#lI-VEVk-6#lI?? S3=S3+V EVK#4I-V EVK-6#4T 

7020 S4=S4+V EVk# 1I»V EVk# 11“V EVk-6# 1I»V EVK-6# IT 

7050 S5=S2»S2?? BN=6^S4-S5?? IF BN=8 THEN CU=CU+1?? RETURN 

7070 M1=<6^S1-S2#S3>/BN?? Tk=M1^68 

7030 H1=<S3»S4-S1^S2>-''BN?? TT=V EVk# 1T»M1+R1 

7110 21=TT?? 22=8228H?? GOSUB 4888?? Zl=Tk?? 22=82R8H?? GOSUB 4888 

7130 MWDE89FF0MT=e.-3^T+3O0?? MWD E89FF2HT*TT«3? ? RETURN 

7200 S1=0?? S2=0?? S3=0?? S4=0?? S5=8 

7210 FOR IS=Vk-7 TO Vk?? Sl=Sl-»-V EIS# 1J#V EIS#4J 

7220 S2=S2+VEIS#1T?? S3*S3>V E1S#4T ? ? S4»S4-**VE1S# IT^V EIS# IT 

7230 NEXT IS 

7240 GOTO 7850 
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3.0 COMPUTER HARDWARE 


The computer is a Texas Instrument TMS990/101M micro-com- 
puter using the TMS9900 16-bit microprocessor and 16K bytes 
(of 8 bits) of RAM. The system presently includes: 

a. TM 990/101M. The 101-M board (slot 1) holds the 
TI9900 microprocessor, a 9901 CRU processor, 4K of 
EPROM (the BASIC interpreter) 2K bytes of RAM and 
various drivers and buffers. The jumper settings 

are El-2, 4-5, 7-8, 9-10, 13-14, 16-17, 18-19, 22-23, 
24-25, 26-27, 28-29, 31-32, 33-34, (these four pairs 
designate type 2716 EPROMs), 39-40, 54-55. Dip switch 
is irrelevant since there is only a single user. The 
two 47 uf capacitors (C8 and C23) debounce the switches. 
Memory on this board maps EPROM from 0000 to IFFFH, and 
RAM from 0F000H to 0FFFFH. The KSR terminal is connected 
by cable to P2 . 

b. TM 990/201-42 (Slot 2), memory expansion board has 
4K RAM and 8K EPROM of which 2K holds the MIAMI program 
and 2K a duplicate. The on-board dip-switch is set to 
l-O-O-l-O-l-l-O to map EPROM onto addresses 2000H to 
5FFF and RAM onto addresses OCOOOH to ODFFFH. Both J1 
and J2 are set to "slow". 

c. TM' 990/302 software development board (slot 3), 

with 4K EPROM holds the development software on EPROM 
and provides 2K RAM. The EPROM is disabled when the 
system is in use. It is used for assembly lainguage 
programming, writing EPROMs and the like. The dip 
switch is set 1-on, 2-off, 3-off to map RAM to ad- 
dresses 0E000H to 0EFFFH. Rocker switch is depressed 
on the side corresponding to the word. Jumper settings 
are E1-E2, E5-E6. The jumper table is standard. The 
cassette recorders are attached to P—2, P-3 holds the 

EPROM "personality" module when "burning" EPROMs. 

d. TM 990/310 input/output expansion board (slot 4), 
provides three identical 9901-driven parallel 

I/O ports (40-pin edge connectors.) Two, P2 and P3, 
are dedicated to driving the panel digital displays 
and one is reserved for future development of inter- 
rupts, and for digital signals out, e.g., turn on de- 
icers, The dip switch is set to On-On-On-On which 
establishes the CRU address of P2 as 0200H, P3 as 280H, 
and P4 as 0300H. P2 delivers the ice thickness number 
to the lefthand PDD, while P3 delivers accretion rate, 
each through 40-wire ribbon cables. 
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MIAMI 












e. Two Analog Devices 1241-R Analog to Digital converters 
are in slots 7 and 8. These boards convert digital data 
to analog voltages or currents for transmission of the 
frequency index to the MIAMI, and signals to external 
displays (scope or plotter) and accept the MIAMI response 
voltage and convert it to a number. The upper board 
provides a gain of 1, and the lower board a gain of 3. 

The input signal is delivered to each in parallel and 
depending on its strength the converted number is read 
from board 1 or 2 (See program.) Each board is jumpered 
to appear as an unoccupied block of 16 bytes of main 
memory of the lOlMi Board A has its block beginning 
at memory address 8ffo, while board B's block is at 
9ffo. (Zeros are strapped on the tables where XI is the 
high order 4 bits and the low order 4 bits are always 
zero.) Each board has 2 DAC outputs, each available as 
both current and voltage and 32 single-ended analog 
inputs. Connector wiring is given in the following table 
followed by a table of measured voltages at various 
points in the transducer. 
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Cable to MIAMI transducer 


Pin 



i (current loop voltage to MIAMI) (gray ribbon) 





Ln 


+ 15 V (green) 
- 33 / (red) 


+15 ✓ (white) 

GND (black, (ribbon) 

+33 / (blue) 

Ai (voltage) from detector amplifier to ADC 
input) (ribbon) 

GND 


Distribution of cable at MIAMI 
Component side, 

+15, -15, A^ (out), i (in) +33, GAD, -33 
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4.0 MIAMI OPERATING INSTRUCTIONS 


1. After installing the transducer, connect the transducer cable, 
oscilloscope cables (top to 'x' , bottom to 'y',) and ground wire 
to scope, plug in MIAMI ^ terminal, scope and connect printer 
cable to printer. 

2. Press latching off /on button on MIAMI power supply. Button 
should show red light. 

3. Press OR on printer, which will type banner message, and 
"ready*' . 

4. Type NEW 0C800H . Banner will print again. 

5 . Type LOAD 4000H . 

6. Type RUN . MIAMI is now operating. Printed matter should be as 

follows: TM996 BftSlC REV D.1.10 

♦RERDY 
NEW 0C300H 

TW990 BRSie REV D.1.10 

♦RERBY 

LDRD 4000H 

RUM 

7. Computer now initializes itself and when ready, sends the 
characteristic resonant sweep to the oscilloscope- Adjust 
scope to place peak in upper right hand corner and left tall 
at left lower corner. During this time, digital displays on 
panel (PDD) will show, left to right, peak amplitude and loca- 
tion of peak, as number of counts. Ignore the decimal points. 

With no protective layer, left should read about 1949 and right, 
about 1875. With 25 mils of plastic, counts will be about 
1423, 1550. Nxombers may jitter in the least significant digit. 
This is due to inherent noise in the system, and not important. 
This completes the test of the system. 

8. When satisfied with the display, type C (CR) . for "cali- 
brate". Displays will read zero. Scope will show a bright 
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peak at top of curve for a moment, then four dots, one in each 
corner of the raster. These dots locate the corners of the ice 
thickness graph which will be generated on the scope as the test 
progresses. Time from 0 to 240 seconds progresses from left to 
right across the scope. Ice thickness, in mils, from 0 to 200 
progresses from bottom to top. Ad:j.ust the scope as required. 

Then erase the display. The dots will return. System is now ready 
and waiting for ice. As ice accumulates, its thickness as a func- 
tion of time will be displayed as a diagonal line on the scope. 

9. When about 2 mils of ice have accumulated on the transducer, 
the beeper will sound and readings of ice thickness will be begun 
by the computer and stored. Begin transit readings at the same 
time. Observer should set transit 20 grads ahead of the ice, call 
out "Mark!" as visible edge of ice passes the transit graticule 
and then move the graticule another 20 grads. At the same time, 
the MIAMI operator should press CR on the terminal. Computer will 
note and store the time at which this occurred for later conver- 
sion to ice thickness as a function of number of marks observed. 

When CR is entered, printer will echo (print) a ".". 

10. The program will run for about 210 seconds, recording two 
measurements each second for the first 30 seconds, then one mea- 
surement per second, until 240 measurements have been made or the 
run is aborted. The run may be aborted by typing "B" any time after 
it has started. The program jumps to the Report Section upon entry 
of "B" or when 240 readings have been made. 

11. The terminal will ask for Run number. Tunnel parameters (vel- 
ocity, temperature, droplet size, water content), comments or ad- 
ditional data, and will print the transit calibration coefficient. 
This number is the number of mils of ice thickness represented by 
20 grads on the transit, and should be predetermined as part of 
the transit calibration procedure. 
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12. The printer will then output a table giving TIME, in sec- 
onds, at which transit reading was noted, D INDEX, a number pro- 
portional to shift of resonant frequency from the calibration 
(zero ice) frequency, TRANSIT ICE thickness as determined from 
elapsed time and calibration coefficient, and MIAMI ICE thickness, 
in mils as computed from D-INDEX. If more than 30 "MARK" entries 
are made, the last entry is repeatedly over-written by its suces- 
sors . 

13. All of the data (up to 240 lines) accumulated during the run 
is now printed. The columns give time of reading, in seconds, 
resonant amplitude, a number proportional to the power generated 
by the MIAMI for given ice thickness at resonance; frequency shift, 
the same number as shown in the previous table list for every 
entry; thickness of ice in mils as computed by MIAMI; and accre- 
tion rate in mils per minute.* Accretion rate is the slope of a 
least square fit to ice thickness over the last 6 readings, times 
60. No calculation of accretion rate is made during the first 5 
seconds. The right hand FDD (accretion rate) does not give nega- 
tive rates (shedding). If the rate is negative, the number ap- 
pearing (which will contain blanks) is meaningless. The printout 
may be terminated by typing a "B" (abort) while the listing is 
being made. The line being printed will be completed, then the 
printer will space four lines and ask for the operator's next 
choice, which can be "A" for the next test, in which case memory 

is zeroed and program returns to Step 11 above; "D" to reprint 
the report in its entirety; or'E" to reprint the report without 
header data. 

14. To terminate the session, simply turn off the computer, the 
printer and oscilloscope. 


* The first row gives in DX the calibration number 
(DX for foo), all the other numbers in this column 
are this number less observed DX, thus being proportional 
to shift in frequency. 
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